Kraus representation for density operator of arbitrary open qubit system 
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We show that the time evolution of density operator of open qubit system can always be described 
in terms of the Kraus representation. A general scheme on how to construct the Kraus operators for 
an open qubit system is proposed, which can be generalized to open higher dimensional quantum 
systems. 
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INTRODUCTION 

It is well-known that for a closed quantum system, its 
time evolution can be described by a unitary operator. 
However, for an open system, the time evolution is not 
necessarily unitary. The evolution of an open system is 
usually described by the Kraus representation Since 
a real physical system is generally entangled with its en- 
vironment, the proper understanding on the nature of 
the Kraus representation for an open s^tem is impor- 
tant and useful 1, 2, 3, H H 0, 0, lE H E3, especially in 
quantum information processing. 

The Kraus representation of an open system is usu- 
ally constructed by considering a larger closed system 
denoted as Sie, comprising of the interested system Si 
and its environment Se- Let pie{t), Pi{t) and Pe{i) be 
the density matrices of Sie , Si and Se respectively, where 
Pi[t) = ive[pte{t)\ and/9e(i) = ivi[pie[t)\, and Pie (0), pi(0) 
and Pe(0) represent the corresponding initial states re- 
spectively at i = 0. As the combined system is a closed 
one, its evolution is unitary. 



p^e{t)^U^e{t)p^e{Q)U^e{t)+, 



(1) 



where Uie{t) is the unitary operator. The interested sys- 
tem, as an open one, then evolves in the following way 



p^{t) = ire{U^e{t)p^e{Q)U,e{t)+}. 



(2) 



If the above equation can be equivalently expressed in 
the form 



p,{t) = ^Af^,(<)p,(0)AV(t)+, 
where Mfj^^{t) satisfy 



(3) 



(4) 



it is said that the evolution of pi{t) has the form of the 
Kraus representation. 

It is obvious that pi{t) always has the Kraus represen- 
tation for arbitrary Uie{t) if pie{0) is factorable p|, i.e. 



Pie{,Q) = Pi(0) (8) /3e(0), which mcans that there is no ini- 
tial correlation between the open system and its environ- 
ment. To show this, we can take Pe{0) = \/PvWe) i^^el 
and let 



(5) 



where |/ie)i We) { Pi = 0, 1, ...,k— 1) are the orthonor- 
mal bases of Se, and k is the dimension of Se, one will 
find that Mf^^{t) defined by Eq. satisfy Eqs. l|2Jl and 

The issue is that whether pi{t) still has the form of 
the Kraus representation when pie(O) is not factorable, 
which means that the initial correlations between Si and 
Se are present. Or in other words, can one always find 
the Kraus representation of an open system for arbitrary 
initial state /Oie(O) and arbit rary unitary operator Uie{t)7 
Recently, some papers 0, S^^JS have contributed to the 
issue. Stelmachovic et al. [g] investigated the role of the 
initial correlations between the open system and its envi- 
ronment and showed that a map based on the reduced dy- 
namics in the presence of initial correlations can't be de- 
scribed by the form of the Kraus representation because 
an additional inhomogcneous part appears. Salgado et 
al. JJ] pointed out that pi{t) still has the Kraus repre- 
sentation even in the presence of any initial correlation if 
the evolution is local, namely Uie{t) = Ui{t)®Ue{t). In a 
very recent paper ^lOj , Hayashi et al. examined the valid- 
ity of the Kraus representation in the presence of initial 
correlations and concluded that the dynamical map for 
an open system reduced from a combined system with an 
arbitrary initial correlation takes the form of the Kraus 
representation if and only if the joint dynamics is locally 
unitary. 

To arrive at the above conclusion, an operator /Ocor (O), 
called the correlation operator, was introduced through 
the definition Pcori'i) = Pie(O) — Pi(0) ® Pe(0). Equation 
^ can then be recast to the following form 

p,(t) = iVe{U^e(t)pM®Pe{Q)U,e{ty} 
+tre{U^e{t)Peor{0)U^e{t) + } 

= Y,M,,,{t)p,{{))M^,{t)+ + 5p,{t), (6) 
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where 



5p^{t) = tre{U^e{t)pcor{Q)U^e{t)+}. 



(7) 



The analysis in Refs. |^ 0, 0| is based on the idea 
that Pi{t) has the Kraus representation if and only if 
5pi{t) = 0. Clearly, pi{t) has the form of Eq. @ if 
5pi{t) = and the Kraus operators are given by Eq. 
(O. However, noticing that Kraus operators are highly 
nonunique, one may start wondering whether pi{t) has 
an alternative form of the Kraus representation even if 
5pi{t) 7^ 0, because there may exist Kraus operators 
Mfj,^{t) such that 

p,{t) - ^AV(i)/O»(0)AV(t)+ + 5p,(t) 

= Y^M,,,{t)p,{0)M^,{t)+, (8) 

and = I- Mf^^ may not be calcu- 

lated from Eq. (|SJ), but they need to have the proper- 
ties of Kraus operators, which ensure the map defined by 
them to be hermitian, trace preserving and positive. 

We consider this problem in the present paper. Our in- 
vestigation focuses on the open qubit system. The paper 
is organized as follows. In Sec. II, an example is pro- 
vided to show that the alternative Kraus representation 
really exists even if 6pi{t) ^ 0. In Sec. Ill, we propose 
a general approach on how to construct Kraus operators 
for an arbitrary open qubit system. We end with some 
discussions in the final section. 



In the model considered, p^e is a 4 x 4 matrix while pi 
and Pe are 2x2 matrices. For simplicity, we take the 
initial state of the combined system as 



/ 



P.e(O) 








V 



(10) 



where tq G (0, 1) is a real parameter. Noting that at tq = 
or 1, Pie(O) is factorable, and the Kraus representation 
certainly exists, we needn't consider these two cases. 

It is easy to obtain the initial reduced density matrices 
of Si and Se as 



Pi{0) = irePie(O) = ^(1 - roa^), 

Pe(0) = iriPie(O) = ^(1 - roCTz), 



and the correlation operator is 



Pcor{0) = ^(1 - ''o) ® CTz- 



(11) 



(12) 



From Eqs. (0), ® and H10() . we get the density matrix 
of the system Si, 



. - _ 1 / 1 + sin^ t — ro cos^ t 
^'^-^^2 1 i(l -Fro) sin i cost 



i{l + To) sint cost 
(1 + To) cos^ t 



(13) 



KRAUS REPRESENTATION WITH NONZERO 

In this section, by providing an example, we show that 
Pi{t) may still have an alternative form of the Kraus 
representation even if 5pi{t) ^ 0. We choose the same 
model as that Ref.^3| has used. That is, we consider a 
combined system composed of two spin- 1/2 subsystems 
with the interaction Hamiltonian Hie — crx^'^{l — (Jz) + 
1 ^{1 + az) , where ax and az are Pauli spin oper- 
ators. In this model, the first qubit plays the role of 
the open system while the second qubit plays the role 
of the environment. The interaction described by the 
Hamiltonian corresponds to the well-known controUed- 
NOT gate|ll3. The unitary evolution operator is given 
by Uieit) = e-'""\ explicitly 





cost 

e-" 

-isint 



-i sint 


cos t 



(9) 



Substituting Eqs. © and (fT^ into Eq. (Tj), one gets 



r f-L\ 1/1 2\ / 2 sin t 



1 sin 2t 

2 sin^ t 



(14) 



We see that 6pi{t) is, in general, non-zero. However, the 
Kraus representation of pi{t) still exists. One can verify 
that the following expressions hold. 



,(t) = ^A/^(t)p,(0)Af^(t)- 



(15) 



(16) 



with 
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Mo(t) 
Miit) 



'\J + ^o)(7't + sin^ t — tq cos^ t) iy' (1 — )(rt - sin^ i + tq cos^ t) 
1^/(1 + ro)(rf - sin^ t + ro cos^ t) \ {I - rt){rt + sin^ t - cos^ i) 



v/2rt(l + ro) 



\/ Tt + SVC? t ~ rQ cos^ i 
i\/rt — sin^ i + t'o cos^ t 



(17) 



where rt = y sin^ ^ + ''o '^'^^^ ^' ^'^o(i) and Mi{t) are the 
Kraus operators. 

The map defined by Eq. (|17|l ensures Pi{t) hermitian, 
trace preserving and positive. The evolution of the sys- 
tem 5*^ obeys Eq. p5|) while the combined system evolves 
under unitary operator Uie{t) given by Eq. ((SJ. This ex- 
ample has showed that even if Spi{t) ^ 0, pi(t) can still 
be written as the form of Kraus representation. 



KRAUS REPRESENTATION FOR ARBITRARY 
DENSITY OPERATOR 

From Eq. lO, we see that the state Pi(t) cannot be 
written in the form of the Kraus representation with the 
Kraus operators defined by Eq. © if 5pi{t) ^ 0. How- 
ever, the example in section II illustrates that there may 
exist an alternative form of the Kraus representation even 
if Spi{t) 0. This encourages us to conjecture that the 
time evolution of the density operator can always have 
the Kraus representation irrespective of the forms of ini- 
tial state and evolution path. In this section, we will 
prove that it is true that pi{t) always can be connected 
with its initial state Pi(0) by Kraus operators. 

Let us begin by considering an arbitrary evolution of 
an open qubit system with arbitrary initial state. The 
most general initial state for an open qubit system can 
be written as 



P^(0) = i(l + ro-a) = ^ 



1 



ro cos f^o 

100 



ro sm Wo e 



rosm^o e"'™ 
1 — ro cos ^0 



and the most general evolution of the system reads as 



p,{t) = -{l + r- a) = 



1 / 1 



r cos f 



r sm f 
1 - 7 



-iff) 



CIS) 



where r = r(t), 9 = 9{t),(j) = </'(i), depending on time 
t, and r(0) = ro, 9{0) = 9q, (/)(0) = 0o- < r < 1, < 
9 < TT, < (p < 2tt. We want to show that there always 
exist the Kraus operators M^{t), such that 



p,(t)=^Af^(<)p,(0)Af^(<)^ 



(20) 



(21) 



where we have used M^{t), instead of M^^{t), to denote 
the Kraus operators. To find the Kraus operators, one 
may write M^(t) as 2 x 2 matrices with undetermined 
elements and one may then directly solves Eqs. (|20|l and 
(|21|l to determine the matrices. However, it is too difh- 
cult to do in that way. As diagonal matrix is, in general, 
easier to handle than non-diagonal ones, we first diag- 
onalize the density matrices Pi{0) and Pi{t) by unitary 
transformations, 

PM = Ulp'mUt. P^it) = U2{t)p[{t)U2{t)+. (22) 

The eigenvalues of Pi(0) and pi{t) make up the entries 
of the diagonalized matrices Pi(0) and p-(t) respectively. 
And their orthogonal vectors make up the columns of the 
unitary matrices Ui and U2 respectively. In this way, the 
diagonalized matrices can be written as 



ro 



cr = - 

' 2 



1 + r 






1 - r 



(23) 



(24) 



where rj, and r' are defined as rj. 



(0, 0, — ro) and 
r' = (0, 0; r) respectively, and the corresponding unitary 
transformation matrices are 



t/2 = 



sm 



(25) 



(26) 



m 

If we can find such operators M'^{t) that satisfy p-(t) = 
A^;WpUO)M;(i)+ and ^^{1)+ M'^{t) = I, then, 
the Kraus representation of pi(t) can be realized by let- 
ting 



M^(t) = U2M'{t)U+. 



(27) 



Since p'^{t) and p^(0) are diagonal, the operators M'^{t) 
are easy to find. There are infinite choices of this kind 
of Kraus operators. Without loss of generality, we may 
choose them as 



Am = 



Mm 



1 







l-r 









(28) 
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Substituting Eqs. if^ . flE^ and (^5)1 into Eq. (|?7j) . we obtain the Kraus operators M^{t), 
I 




00 

2 



sm 2 cos -fe^ 



COS I COS ^e-*"^" - 



J:-^ sin I sin ^e"'"^ 

l+io 2 2 



COS I COS 2f e*"Po 

1+ro 2 2 



COS ^ sm ^ 

2 2 



cos I cos fe'f" 
sin f cos |te*("^+'^°) 



sm 5 sm ^e'"^ 



(29) 



Mo{t) and Afi(i) satisfy Eqs. ^ and (ETJ- The Kraus 
representation given by Mo(t) and AIi{t) in expression 
(|29|l does ensure Pi{t) be hermitian, trace preserving and 
positive. So, no matter what the forms of Uie{t) and 
Pie(O) are, there always exist the Kraus operators con- 
necting piit) with Pi(0). For any given Pie{0) and Uieit), 
the Kraus operators (i) can be calculated by diagonal- 
izing the reduced matrices ts:e.{Uie{t)pie{^)Uie{t)^} and 
tre/3ie(0). One general expression of the Kraus operators 
is given by Eq. H29() , with which the Kraus representation 
of pi{t) is obtained by Eq. 

So far, we have proved that the time evolution of an 
density operator of open qubit system always has the 
Kraus representation. At the same time, we have put 
forward a general approach for constructing the Kraus 
operators for arbitrary evolution. From physical point of 
view, the above process of finding the Kraus operators 
means that we first align the Bloch vectors r and Yq in 
Bloch sphere along the z axis by using Ui and U2 respec- 
tively, find the Kraus representation of r' with rj,, and 
then revers r' and Tq back to r and Tq, to obtain the 
Kraus representation of r with Yq. The model in sec. II 
is just an example of applying this approach to solve the 
Kraus representation. In fact, expression (|17(l is calcu- 
lated in this way. 

DISCUSSIONS 

We have shown that the time evolution of the density 
operator of an open qubit system can always have the 
Kraus representation. A scheme on how to construct the 
Kraus representation is proposed. One general expression 
of the Kraus representation for an arbitrary evolution is 
provided by Eqs. (|2U and Since the expres- 

sions of the Kraus operators are not unique, the form 
given by Eq. H29|l is only one kind of them. The equiv- 
alent expressions of the Kraus operators can be written 
down as M^{t) — Ylv ^'^i^)^i^>'^ where V^^ are the ele- 
ments of an arbitrary unitary matrix. 

Refs. 1^ 0, 0| have investigated the possibility of 
the Kraus representation for an open system with ini- 
tial correlations between the system and its environment 



and some important conclusions have been derived. As 
a supplement, the present paper studies the existence 
of an operator-sum representation for an arbitrary given 
evolution of density operator. Our result shows that an 
arbitrary evolution of the state can always be written as 
the form of the Kraus representation. The Kraus opera- 
tors can be calculated by Eq. ||SJ| if Spi{t) = 0. However, 
they cannot be expressed explicitly in the form of Eq. 
Q if 5pi{t) 7^ 0. For this latter case, they can still be 
obtained by the approach described in the current paper. 
Moreover, M^{t) are generally dependent on the initial 
state and there does not exist a universal form of Kraus 
operators for all different initial states. 

This approach can be generalized to higher dimen- 
sional quantum systems. The procedure for higher di- 
mensional systems is similar to the qubit case but may 
be more complicated. In fact, the density matrix Pi{t) 
with the parameter t and Pi(0) can always be diagonal- 
ized by unitary transformations Ui and U2 respectively, 
regardless of the dimensions of the density matrices. It is 
easy to find the Kraus operators M'^{t) of diagonal den- 
sity matrix p[{t) with diagonal initial density matrix p[{Q) 
though it is difficult to find the Kraus representation of 
an arbitrary density matrix with arbitrary initial condi- 
tions. By Eq. (EUl, using L/i, U2 and M'^[t), the Kraus 
representation of pi{t) is obtained. Certainly, as the di- 
mensions of the density matrices become larger, solving 
for the Kraus operator may become more formidable. 
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